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Abstract 

We present Monte-Carlo simulations for the focusing design of a novel cold-neutron triple-axis spectrometer to be 
installed at the end position of the cold guide NL-1 of the research reactor FRM-II in Munich, Germany. Our simula- 
tions are of general relevance for the design of triple-axis spectrometers at end positions of neutron guides. Using the 
McStas program code we performed ray trajectories to compare parabolic and elliptic focusing concepts. In addition 
the design of the monochromator was optimized concerning crystal size and mosaic spread. The parabolic focusing 
concept is superior to the elliptic alternative in view of the neutron intensity distribution as a function of energy and 
divergence. In particular, the elliptical configuration leads to an inhomogeneous divergence distribution. 
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1. Introduction 

Neutron focusing techniques have become rather im- 
portant in triple-axis spectroscopy because there is a 
high demand for measurements which require an in- 
tense, focused neutron beam on small samples. Dou- 
ble focusing at the monochromator side combined with 
a horizontally focusing analyzer were introduced by 
Biihrer [1] and Pintschovius H in the 1980s and be- 
came very successful. The double focusing monochro- 
mator allowed for inelastic neutron studies of the 
phonon and magnon excitations with the first avail- 
able small crystals of high-temperature cuprate super- 
conductors [3] as well as for the routine measurement 
of the phonon dispersion under pressure up to 10 GPa 
J3-01- Especially the ongoing interest in experiments 
using high-pressure cells with very small samples re- 
quired further experimental progress which could be 
achieved by the use of focussing neutron optics first in- 
troduced in front of the sample by Goncharenkov et al. 
ii. Using elliptical neutron optics it became possible 
to study samples as small as 0.1 mm 3 by neutron scat- 
tering techniques i?- IT], which is of particular impor- 
tance in the case of novel materials like the new, ex- 
citing class of high-temperature superconducting iron- 
pnictide compounds Il2l - ll5ll . as it it is frequently dif- 
ficult to obtain large samples. Focused neutron beams 
are also of general importance for inelastic neutron scat- 



tering experiments since the gain of intensity results 
in a decrease of measuring time, in an increase of the 
measurement-statistics, and often also in a better signal 
to noise ratio. 

Due to the ongoing progress in neutron optical de- 
velopments it is now possible to use supermirror guides 
with different shape or tapering designs and rather large 
angles of reflection of the supermirror coatings 11161- 
l20ll Together with the virtual source concept |B l2lll 
and large focusing monochromator arrays [|22h26|| the 
intensity at the sample position can be enhanced signif- 
icantly. 



Recent simulations [ 27] indicate a flux gain of the or- 
der of five which could be obtained by using an elliptic 
focussing guide instead of a conventional guide with an 
m-value of two (m is a measure of the reflection angle 
relative to the angle of reflection of Ni). 

Here, we compare elliptic, parabolic and conven- 
tional (straight) guide concepts for a cold-neutron triple - 
axis spectrometer using the free McStas code l28T - [30ll . 
Our results show that for an instrument at the end of 
a cold neutron guide, the parabolic concept is clearly 
superior to the elliptic concept concerning the intensity- 
divergence distribution. 
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Figure 1: Configuration of the cold-neutron spectrometer. First, a 
parabolic guide is installed permanently. The following two parabolic 
guides can be removed and replaced by other guide elements. 



2. The optimized focusing concept for a cold- 
neutron triple-axis spectrometer 

Our calculations were made to design the concept of 
a new cold-neutron triple-axis spectrometer which will 
be installed at the end of a curved neutron guide (NL1) 
at the research reactor FRM-II in Garching. Therefore 
we have used the specific dimensions and coatings as 
well as the intensity distribution of this neutron guide, 
but our main results are of broad relevance for any cold- 
neutron triple-axis spectrometer placed at a neutron- 
guide end-position. A gap of 400 mm is required for 
the monochromator of the upstream instrument. The 
width and the height of NL-1 (coating m — 2) amount to 
60 mm and 120 mm, respectively. In Ref. 113111 a fully el- 
liptic guide concept has been presented and a distinctly 
larger entrance width (70 mm) of the guide after the po- 
sition of the upstream instrument has been proposed in 
order to overcome the neutron loss caused by the gap 
where the monochromator is hosted. 

Fig. Q] shows the overall design of the calculated 
parabolic focussing concept of the triple-axis spectrom- 
eter. The parabolic guide has a total length I = I1+I2 + I3 
of 6.7 m and is divided into three sections h, h and I3. 
The first section l\ has a length of 5.15 m and is in- 
stalled permanently. The other two sections I2 and /3 
with 1 .4 m and 0. 15 m length can be moved transversely 
and replaced by other (straight) guide elements of iden- 
tical length. The largest section l\ also hosts the polar- 
izing cavities. Its design will be subject of a forthcom- 
ing publication [32]. Several focal distances have been 
studied in great detail ranging from 0. 1 m to 0.6 m. Fi- 
nally, a focal length / = 0.3 m has been chosen as an 
optimal value for highest resolution and intensity. After 
the parabolic guide a velocity selector will be installed 



in order to suppress the higher harmonics and to reduce 
the background of the instrument. Additionally, sev- 
eral diaphragms and flight tubes shall be installed be- 
tween the parabolic guide and the monochromator drum 
for similar purposes (e.g. background reduction). An- 
other diaphragm at the focal point may also be used as 
a virtual source if even higher energy resolution is de- 
sired. In the monochromator drum a double focusing 
monochromator is installed at a distance dpM away from 
the focal point; compare Fig. Q] The size of the highly 
oriented pyrolythic graphite (HOPG) crystals has been 
chosen as in Ref. 13111 . i.e. with 20 mm x 20 mm size, 
and the mosaic spread is assumed to be 0.5°. All con- 
figurations have been compared based on these identi- 
cal conditions. Finally, we will also discuss a differ- 
ent monochromator design which is optimized for the 
short distances dp u and cIms = 1.2 m with d/ws being the 
distance between monochromator and sample; compare 
Fig. Q] Both, the intensity at the sample position and 
the peak profile as a function of horizontal divergence 
and wavelength could be altered and, hence, designed to 
our needs by variation of the distance dpM — dMs ~ 

Ad. 

At the sample position a neutron counter of 20 mm x 
20 mm size has been used in order to measure the neu- 
tron flux at the sample position. 



3. Comparison of different focusing concepts 

Basically, there are six different focusing configura- 
tions which will be discussed in this work: 

• parab. #1 parabolic guide {l\ to Z3 parabolic) 

• parab. #2 parabolic guide (^3 is straight guide) 

• parab. #3 parabolic guide (I2, h are straight) 

• elliptic #1 elliptic guide (all noses; / = 0.3 m) 

• elliptic #2 elliptic guide (no noses; / = 1.2 m) 

• straight standard straight guide 

The first configuration, parab. #1 , is a fully parabolic 
guide with / = 0.3 m. In the second parabolic config- 
uration, parab. #2, the last nose with 0.15 m length is 
replaced by a standard straight guide. In parab. #3 also 
the much longer second nose I2 is replaced by a straight 
guide changing the properties of this concept distinctly. 
For the sake of simplicity, the third parabolic concept 
will be discussed later in section 4. All parabolic config- 
urations start with 60 mm width of the entrance window 
after the gap for the upstream instrument. In contrast, 
for the two elliptic configurations an increased entrance 
width of 70 mm has been applied since this widening 
of the guide was recommended in Ref. [31]. Regarding 
the two elliptic setups, all elliptic noses were used in the 
first elliptic configuration, elliptic #1, with / = 0.3 m, 



2 




5.08 





guide 




ellip. #1 




ellip. #2 


— • — 


parab. #2 




parab. #1 



5.24 



5.4 




Parameter 


parab. #1/ 
parab. #2/ 
parab. #3 


elliptic #1/ 
elliptic #2 


straight 


Entrance width 


60 mm 


70 mm 


60 mm 


Exit width 


12.4 mm/ 


27.8 mm/ 


60 mm 




15.2 mm / 


52.8 mm 






21.6 mm 






Entrance height 


120 mm 


120 mm 


1 20 mm 


Exit height 


120 mm 


120 mm 


120 mm 


focusing guide length 


6.70 m/ 


6.70 m/ 


m 


(length of parabolic 


6.55 m/ 


5.80 m 




or elliptic guides) 


5.15 m 






straight guide length 


ml 


0m/ 


6.70 m 


( of guides following 


0.15 m / 


m 




the focusing guides) 


1.55 m 






Coating (hor.) 


av. 3.5 


4 


4 


Coating (vert.) 


2 


2 


2 


focal distance / t 


0.3 m/ 


0.3 m/ 


oo 




0.45 m/ 


1.2 m 






1.85 m 






focal distance /' i 


oo 


2.5 m 




Distance gap-F ft 


7.0 m 


7.0 m 


7.0 m 


Mono, shift Ad 


-0.3 m 


-0.1 m 


-0.1 m 


dins 


1.2 m 


1.2 m 


1.2 m 


Rhor 


1.4/1.4/1.6 


1.3 





Figure 2: Neutron intensity at the sample position as a function of 
wavelength for the five different guide concepts (see text). 



Table 1: Guide parameters for the five different focusing concepts, f: 
/ is the distance from the end of the last parabolic or elliptic com- 
ponent to the focal point. Only horizontal focusing has been applied 
for all models. %: f is the focal distance of the entrance of the ellip- 
tic guide. J: 'gap-F' is the distance between the gap of the upstream 
instrument and the focal point F. Rhor - Horizontal monochromator 
curvature prefactor (see text). 



and all noses were removed in the second elliptic con- 
figuration, elliptic #2, yielding / = 1.2 m. Finally, the 
configuration straight is a standard straight guide with 
length I. The distance between gap and focal point is 
identical for all concepts compared in this work. All 
guides have supermirror-coatings with m = 2 and 4 for 
the top/bottom and the sides, respectively (see Tab. [TJ. 
However, for the final parabolic configuration an addi- 
tional intensity gain of about 10% could be achieved by 
an optimized distribution of the m-value of the super- 
mirror coating along the parabolic guide with smaller 
values at the beginning but larger values at the end. 
Thus, the guide starts with an m-value of 2.5 which con- 
tinuously increases up to a value of 3.5 at the end of sec- 
tion l\, i.e. after 5. 15 m length. Only the two short noses 
h and /3 possess coatings with higher m-values which 
continuously increase from 4 at the beginning of I2 to 6 
at the end of nose h. Hence, on average, the m-value of 
the supermirror coating for the parabolic configurations 
amounts to 3.50. As mentioned before, compared with 
a uniform coating with m = 4 the average intensity gain 
is of the order of 10%; it is smaller for lower neutron 
energies and only slightly higher for very high neutron 
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Figure 3: (a) Intensity (amplitude) and (c,c') energy resolution 
(FWHM) for the five different models parab. #1 (violet circles), 
parab. #2 (yellow-orange circles), elliptic #1 (green triangles), el- 
liptic #2 (red triangles) and straight (blue squares) as a function of 
kj. (b) Ratio of the intensities of parab. #2 and elliptic #1 (closed 
symbols) or straight guide (open symbols), (d) Ratio of the energy 
resolution (FWHM) of parab. #1 and elliptic #1 (closed symbols) or 
the straight guide (open symbols). 



energies. 

Regarding the vertical monochromator curvature the 
standard monochromator curvature (as expected for fo- 
cusing from focal point to the sample position by the 
Rowland condition) has been chosen for all focusing 
guides. For the horizontal monochromator curvature 
a different value has been chosen which has been op- 
timized by comparison of intensity and resolution for 
different horizontal monochromator curvatures in 10% 
steps from 100% of the nominally optimal curvature 
(Rowland condition) up to 200% of this nominal value. 
Finally, a slightly relaxed curvature around 140% of the 
nominal value turned out to be ideal regarding the inten- 
sity. Thus, the standard horizontal curvature has been 
multiplied by a factor Rhor which is given in Tab.Q] 
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Figure 4: Neutron intensity directly after each guide as a function of 
energy and horizontal divergence for the parab. #2, elliptic #1 and 
Straight concept. The intensity / is plotted logarithmically ranging 
from / < 10 5 (blue) to / = 10 9 (dark red). 



For the elliptic models, the guide parameters from 
1 3lj, |33J have been taken without any further optimiza- 
tion. (For the other models, all parameters have been 
always optimized for highest intensity at the sample po- 
sition!) In the elliptic case the total length of the elliptic 
guide is defined by the focal distance /. For the ellip- 
tic #1 concept it is identical to parab. #1 concept and 
for elliptic #2 it is 0.9 m shorter. For the elliptic cases 
/' is the other side focal point at the entrance window; 
compare also IB3I1 . 



3.1. Neutron intensity and energy resolution 

Fig. [2] shows the neutron intensity calculated at the 
sample position as a function of wavelength for the five 
different configurations and for different values of the 
incident neutron wave vector k, ranging from 1.2 A" 1 to 
3.2 A" 1 . Additionally, the intensities for elliptic #1 and 
elliptic #2 have been calculated but without the guide 
widening to 70 mm proposed in reference lUlll . i.e. with 
the standard entrance window of 60 mm width only 
(green and red lines in Fig. 0. The widened guides 
are not useful. On the one hand more neutrons are col- 
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Figure 5 : Intensity-ellipsoids at the sample position as a function of A and hor. div. for the different concepts (see text) at different values of kj . The 
intensity is plotted from 0.0 (blue) to 3 ■ 10 6 (dark red) in arbitrary units, (parab. #3 and parab. #2b will be discussed in the sections 3 and 6.) 



lected by the wide -spread entrance window of such a 
widened guide, but it becomes more difficult to focus 
such a beam on a small 2 cm x 2 cm spot at the sample 
position once it is initially widened. 

As can be seen from Fig. [2] the energy resolution 
of the two parabolic configurations is distinctly better 
than that of the elliptic configurations. Furthermore, 
the neutron intensity of the guides with parabolic taper- 
ing is distinctly larger, especially for larger values of k,. 
In Fig. [3] a quantitative evaluation and comparison of 
parabolic and elliptic concepts is presented showing the 
intensity and energy resolution ratios of both concepts 
(closed stars). The parabolic concept performs already 
better for small incident neutron energies; the energy 
resolution, i.e. the peak width (FWHM), of the com- 
peting high-resolution elliptic setup elliptic #1 is about 
one third larger (worse) than that of the parabolic setup 
parab. #1 . For higher neutron energies, the difference 
is even more dramatic since the energy resolution, i.e. 
the peak width (FWHM), of elliptic #1 has become even 
two thirds larger (worse) than that of the parabolic setup 
at a value of kj equal to 3.2 A -1 (~21 meV). Further- 
more, the neutron intensity at this value of kj has be- 
come almost two times larger than that of the elliptic 



alternative; see Fig. [3] 

In Fig.[3](b,d) a comparison of the parabolic concept 
with the simple straight guide is shown (open stars). 
Whereas the parabolic concept is always distinctly bet- 
ter than a straight guide, the elliptic alternative 1131 1 
loses most of its advantages at higher neutron energies; 
compare also Fig. [2] 



3.2. The transverse Q-resolution 

In Fig. |4] the neutron intensity is plotted as a func- 
tion of horizontal divergence and neutron energy. For 
the standard straight guide the horizontal divergence of 
the neutrons is rather small compared to the fully focus- 
ing concepts. In the other extreme, the parabolic guide 
parab. #2 exhibits a broad divergence. From these re- 
sults, the straight guide with the worst energy resolution 
has the best transverse Q-resolution and the parabolic 
guide with the best energy resolution has the worst 
transverse Q-resolution. elliptic #1 looks like a super- 
position of straight and parabolic concepts. Thus, the 
intensity distribution of the elliptic guide is not as ho- 
mogenous as for the parabolic or for the straight guide. 
Hence, with elliptic focusing one might expect three 
major peaks for a transverse scan in a neutron scattering 
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Figure 6: Schematic focusing design of (a) fully elliptic and (b) fully parabolic guides in the horizontal plane. (a',b') The intensity distributions of 
neutrons at the entrance window are shown as a function of divergence and wavelength together with a schematic drawing of the entrance window 
to the right (vertical magenta bars) and the neutron trajectories (black lines) which belong to neutrons that will be focused into the focal point of 
the corresponding focusing guide. The magenta lines connecting the schematic drawing of the entrance window with the contour maps indicate 
from which divergence regime in the contour maps these neutrons originate depending on their position at the entrance window. 



experiment. In order to elucidate this question, we have 
also calculated the corresponding neutron intensities at 
the sample position after scattering by the monochro- 
mator. 

In Fig.|5]the neutron scattering intensity at the sample 
position is shown as a function of neutron wavelength 
and horizontal divergence for the different models men- 
tioned above. The calculated intensity maps are directly 
related to the resolution ellipsoid of the spectrometer. 
As can be seen in Fig. [5] both elliptic concepts always 
exhibit a multiple peak structure. Therefore, also a per- 
fect single crystal may exhibit multiple peaks in a trans- 
verse scan. In contrast, the multiple peak structure is 
absent in the parabolic concept for medium and high 
neutron energies. Only at rather low neutron energies 
a similar multiple peak structure may appear. In the 
parabolic arrangement the multiple peak structure, how- 
ever, can be suppressed when replacing the focussing 
noses by straight guides (parab. #3), see Fig. [5] 

3.3. The qualitative understanding of parabolic versus 
elliptic configurations 
The better performance of the parabolic versus the el- 
liptic focusing concept for an instrument at the end of a 
neutron guide can be understood by the following rea- 
son: the parabolic guide focuses neutrons with small 
divergence from the whole entrance window into its fo- 
cal point, whereas the elliptic guide focuses neutrons 
with higher divergence from outer parts of the entrance 
window into its focal point. This can be seen in the 
schematic drawings of Fig. |6](a,b). At larger distances 



from the reactor, the neutron intensity provided by the 
guide is dependent on the neutron divergence. As can 
be seen in Fig.|6](a',b') the intensity is strongly reduced 
for neutrons with higher divergence due to the reflection 
losses in the preceding neutron guide. This reduction of 
intensity for neutrons with higher divergence becomes 
especially significant for higher neutron energies. Since 
the elliptic concept also depends on focusing neutrons 
out of these high-divergence regimes into its focal point, 
it has a clear disadvantage with respect to the intensity, 
especially for higher neutron energies. In contrast, the 
parabolic concept basically focuses neutrons with low 
divergence (from the regime with highest flux in the 
intensity-divergence distribution) into its focal point. 

4. Modification for High-Q resoiution 

With a basic modification of the parabolic concept 
one may enhance the Q-resolution. As can be seen in 
Fig. Q] it is possible to remove both parabolic noses of 
the calculated design (fa an d '3) an d to exchange them 
by straight guides. We will denote this alternative setup 
by parab. #3. This parabolic concept is almost equiv- 
alent to the elliptic #2 concept. Like for the other con- 
cepts we have optimized the horizontal monochromator 
curvature for neutron flux at the sample position. 

In Fig. [5] the intensity ellipsoids of the different fo- 
cusing concepts parab. #1, parab. #2, elliptic #1, ellip- 
tic #2, straight and parab. #3 are shown for the most 
relevant values of kj. As mentioned before, the elliptic 
concepts exhibit a multi-peak structure for all incident 



6 



neutron energies. The parabolic concepts do not show 
such a multiple peak structure except at low energy. But 
in contrast to the elliptic concept, the parabolic concept 
loses this multiple peak structure when all noses are re- 
moved, i.e. in the parab. #3-design as can be seen in 
Fig. [5] Regarding the intensity at the sample position, 
this concept is comparable to the elliptic concept in the 
lower or medium energy regime but still clearly outper- 
forms the elliptic alternative in the high energy regime. 
This configuration will also have a better Q-resolution 
as discussed in the next section. 

5. Simulation of measurements 

In order to study the different focusing concepts un- 
der measurement conditions, we have also simulated 
several types of scans: 1 ", energy scans for an incoher- 
ent scatterer, 2 nd , elastic scans across a Bragg peak and, 
3 , scans for measuring an acoustic phonon dispersion 
will be presented. 

The double-focusing analyzer was composed of 1 1 
rows and 21 columns of 10 mm x 20 mm large HOPG 
crystals with 0.4° mosaic spread. The sample-analyzer 
and analyzer-detector distances are rather small, i.e. 
0.8 m each, since the spectrometer is intended to have 
a rather compact design. The width of the detector was 
set to 2 cm and all scans have been performed in the 
constant-A:/ mode of operation. Furthermore, no colli- 
mation was used at any point corresponding to an all- 
open configuration. 

5.1. Energy scan for an incoherent scatterer 

First, we have calculated the intensity and energy res- 
olution for energy-scans across the elastic line of a sim- 
ple cylindrical incoherent scatterer of 20 mm height and 
5 mm radius. Each point in reciprocal space was cho- 
sen such that Q was equal to kj and a double-focusing 
analyzer geometry has been chosen. The resulting in- 
tensities and peak widths (FWHM) are shown in Fig. [7j 
As can be seen, the intensity of parab. #2 is always 
larger than the intensity either of the standard straight 
guide or of elliptic #1 . At the same time the peak width 
of parab. #2 is smaller than that for the other concepts. 
Compared to the standard straight guide, the energy res- 
olution is nearly a factor 2 better for parab. #2. It is even 
around 20% better compared to the high-resolution el- 
liptic concept elliptic #1 . Especially at very high energy 
around 30 meV, the resolution of the parabolic concepts 
becomes almost 30% better and the intensity 90% larger 
when compared to the elliptic concept. 

Comparing elliptic and straight guides, the elliptic 
guide always has a clearly better energy resolution. 



Even the intensity is slightly larger for the elliptic con- 
cept at lower neutron energies. However, the intensity 
performance changes at higher neutron energies where 
the straight guide exhibits a better performance. 

5.2. Elastic scans 

Next, we have simulated <y-scans across a (100) 
Bragg reflection of a cubic single-crystal with 5 A lat- 
tice parameter and 50' mosaic spread using a flat ana- 
lyzer geometry. But since it is unfavorable to work with 
a flat monochromator, if the preceeding parabolic or el- 
liptic guide is strongly focusing, we did not change the 
monochromator curvature. The resulting intensities of 
such w-scans are shown in Fig. [8] for two different in- 
cident neutron energies. The highly focusing concept 
parab. #2 exhibits the lowest performance. But also the 
highly focusing alternative elliptic #1 is not competitive 
with a straight guide. Interestingly, again, the straight 
guide exhibits the best intensity-performance for such 
elastic scans across a Bragg peak. However by chang- 
ing towards the parab. #3-configuration even higher in- 
tensities than for the straight guide can be obtained; see 
Fig. [8] Whereas the two very similar fully and almost 
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Figure 7: (a) Intensity (amplitude) and (c) peak widths (FWHM) 
of simulated energy-scans across the elastic line of an incoherent 
scatterer; orange circles: parab. #2, green triangles: elliptic #1 , 
blue squares: straight. In (b) and (d) the corresponding ratios are 
shown; green/orange triangles: ratio of parab. #2 versus elliptic #1 , 
blue/orange squares: ratio of parab. #2 versus straight. 
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fully parabolic configurations parab. #1 and parab. #2 
are optimized for highest energy resolution, the 'semi- 
parabolic' configuration parab. #3 is optimized for Q- 
resolution. 

5.3. Phonon measurements 

Finally, we have simulated inelastic scans across 
transverse and longitudinal acoustic branches of the 
phonon dispersion of a fee single crystal using the 
Phononsimple-component provided by McStas which 
is based on the expressions of cross sections provided by 
Squires 13411 . The model sample has a cylindrical shape 
with 2 cm diameter and 2 cm height and the crystal lat- 
tice has 5 A lattice parameters. A low sound velocity 
corresponding to a slope of the phonon dispersion which 
amounts to 8 meV/A -1 (at 300 K) has been used. The 
sample mosaic is not included in the Phonon simple- 
component, i.e. r\ — 0. For these simulations a double- 
focusing analyzer geometry has been applied. 

The simulated intensities of constant-Q-scans are 
shown in Fig. [9] The intensities have been calculated 
at the detector position. In order to avoid any focus- 
ing effects of the resolution ellipsoid with the phonon 
dispersion (focusing side and defocusing side), we have 
studied first the phonon dispersion at the zone bound- 
ary because its dispersion is flat, i.e. we have sim- 
ulated constant-Q-scans at (1 1 0) in reciprocal space 
with varying energies around ~6.4 meV; compare also 
Fig. [10] (a). We will denote this type of scan by 'scan 
A' and for the sake of simplicity we will only show the 
results for the three major concepts with either (nearly) 
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Figure 8: Simulation of w-scans through (1 0) at 2 meV with (a) 
kf = 2.0 A -1 and (b) kf = 1.5 A -1 respectively. Orange circles: 
parab. #2, green triangles: elliptic #1, red triangles: elliptic #2, blue 
squares: straight, gray circles: parab. #3. 



full or no focusing properties. In Fig. |9](a-f) the result- 
ing energy-scans are shown for different values of the 
final wave vector kf. For a more quantitative analysis 
the peak intensities measured at (Q,ha>) and the corre- 
sponding phonon peak widths are shown in Fig.[9](g,i). 
Finally, in Fig.|9](h,j) the ratios of these phonon widths 
and intensities for the parabolic versus the two other 
models are shown. Obviously, the parabolic concept has 
the highest energy resolution together with the highest 
intensities among the different focusing concepts. The 
energy resolution is 30% better than for the elliptic con- 
cept with at least 50% more intensity. Compared to the 
standard straight concept the energy resolution is more 
than 70% better, having even 100% more intensity in 
the phonon peak. 

Comparing the elliptic and straight concepts, the el- 
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Figure 9: (a-f) Results of constant-Q-scans of type 'scan A' at 
(1 1 0) for parab. #2 (orange circles), elliptic #1 (green triangles) and 
Straight (blue squares) for different kf. (g-j) The phonon peak widths 
and intensities as a function of kf together with the ratios of the dif- 
ferent concepts; green/orange triangles: ratios of parab. #2 versus 
elliptic #1 concept, blue/orange squares: ratios of parab. #2 versus 
straight concept. 
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liptic alternative has larger intensity for medium and 
lower values of kf but loses this advantage of intensity 
versus the straight guide at higher kf. However, the el- 
liptic concept can still keep its advantage in the energy 
resolution compared to the straight guide even at higher 
values of kf. 

In order to study also the focusing effects of the res- 
olution ellipsoid with the phonon dispersion, we have 
simulated constant-energy scans at low energies across 
the zone center for the different basic concepts. The 
scan directions of two constant-E-scans through (2 0) 
are indicated in Fig. [I0](a,b), one of them having a lon- 
gitudinal scan-direction ('scan B') and the other hav- 
ing a transverse scan-direction ('scan C), thus being 
sensitive to phonons with longitudinal and transverse 
polarization, respectively. Both constant-E-scans have 
been calculated with 2 meV energy transfer and in the 
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with kf = 1.5 A -1 , (d) Transverse scan through (200) at 2 meV with 
kr = 1.5 A -1 . Orange circles: parab. #2 concept, green triangles: 
elliptic #1 concept, blue squares: standard straight guide. 



constant-^/ operation mode with kf = 1.5 A~ . As can 
be seen in Fig.[10](c,d) there is no significant difference 
between elliptic (elliptic #1) and parabolic (parab. #2) 
concepts but interestingly the standard straight guide is 
clearly superior to the focusing guides, either elliptic 
or parabolic. This can be attributed to the higher Q- 
resolution of the straight guide. 

Complementary constant-Q-scans at (2.2 0) which 
are not shown here, corroborate these findings - the in- 
tensity of the straight concept is clearly superior to the 
two focusing concepts. However, the energy resolu- 
tion (FWHM) of the different concepts remains almost 
identical and amounts to 0.60(6) meV, 0.59(7) meV,and 

0. 61(6) meV for the parab. #2, elliptic #1and straight 
concepts. Thus, in constant-Q scans measured at a steep 
dispersion, the parab. #2 concept exhibits less intensity 
than the standard straight guide but not a worsened reso- 
lution. Simulations of scans of type 'scan-B' and 'scan- 
C within the parab. #3 concept show no significant im- 
provement, see Fig.[l0](c,d). However, a larger sample 
mosaic will be beneficial for the parabolic setup versus 
the other concepts since larger parts of the divergence- 
band will interact with the sample (remember that zero 
sample mosaic was used in these phonon calculations). 

6. Optimization of the Monochromator 

Because the distance cIms between sample and 
monochromator is rather small for the intended com- 
pact design and for focusing from point to point (with 
Ad ~ 0.3 m), a reduction of the HOPG monochroma- 
tor crystal size seems useful in order to approach a truly 
curved surface since the curvatures become rather large. 
Furthermore, an increase of the mosaic spread could be 
beneficial for these purposes and also for transporting a 
larger part of the divergence-bandwidth provided by the 
cold guide. However, the gap between the monochro- 
mator crystals as well as the required higher precision 
of the monochromator mechanics limit the reduction of 
the crystal size. Enhancing the monochromator crystal 
mosaic might result in a coarser resolution in measure- 
ments with a flat monochromator setting without any 
collimation. Hence, only a small crystal size reduction 
and a small increase of the crystal mosaic seems to be 
useful. After simulating several different setups with 
different crystal sizes (ranging from 10 mm x 10 mm 
to 20 mm x 20 mm) and with different crystal mosaic 
(ranging from 0.4° to 1.0°), we decided to choose a 
monochromator array with 14 mm x 14 mm crystal size, 

1. e. a total reduction of crystal size (area) by ~50%. 
Even smaller crystal size would yield higher simulated 
intensities at the sample position but might result in 
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Figure 11: Calculated neutron scattering intensity at the sample position for different concepts: parab. #2 with optimized monochromator, i.e. 
parab. #2b (black circles) versus parab. #2 (yellow circles), parab. #3 (grey circles), elliptic #1 (green triangles), elliptic #2 (red triangles) and 
Straight (blue squares) with the standard monochromator setup. 



mechanical problems of the monochromator construc- 
tion. Furthermore, a mosaic spread enhanced by 70% 
was chosen as the optimum for the intensity at the sam- 
ple position. We denote this configuration of optimized 
monochromator in combination with the parab. #2- 
concept as the 'parab. #2b' -concept. The combination 
of parab. #3 with this optimized monochromator is de- 
noted as parab. #3b. In Fig.[TT]the results of the simula- 
tions for the high-E-resolution configuration parab. #2b 
(solid black circles) are shown and compared with other 
configurations all having the standard monochromator 
design according to Ref. 13 ill . For all calculated neu- 
tron energies the gain of total neutron-flux at the sam- 
ple position is of the order of one third higher than the 
high flux of the 'conventional' parab. #2 configuration. 
Furthermore, the multiple peak structure transverse to 
Q vanishes since the intensity is smeared out due to the 
larger mosaic spread (see Fig. |5). 

We have also simulated energy scans across the elas- 
tic line of an incoherent scatterer under the same con- 
ditions as described in the preceding section. The re- 
sulting intensities are shown in Fig.Q~2](a) and the ratios 
of the high-E-resolution parab. #2b and the other dif- 
ferent high-E-resolution focusing concepts as well as 



the standard straight guide are shown in Fig. [12] (b). 
Thus, by choosing an optimized monochromator de- 
sign (black solid points), a distinct increase of inten- 
sity can be expected in inelastic measurements in gen- 
eral. Fig. [12] (c,d) show the high performance of 
parab. #2b (black solid points) for measuring flat dis- 
persions. However, Fig. [12] (e) shows that the standard 
straight guide is still superior for measuring steep dis- 
persions. 

Furthermore, we were also curious to see whether any 
performance collapse appears in the parab. #2b con- 
figuration in elastic scans since the crystal mosaic has 
been enhanced. Therefore, w-scans have been simu- 
lated under the same conditions as described in the pre- 
ceding section. In spite of using a flat analyzer geome- 
try there is no drop but even a small increase of inten- 
sity for parab. #3b (solid black stars) as can be seen in 
Fig. [T2](f). 

Hence, the monochromator can be optimized for the 
use in double focusing geometry by choosing crystal 
size and crystal mosaic properly depending on the oc- 
curring monochromator curvatures which are larger for 
a fully focusing guide than for a non-focusing straight 
guide. Since elastic scans usually do not suffer from a 
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Figure 12: (a) Intensity (amplitude) of simulated energy-scans across 
the elastic line of an incoherent scatterer. In (b) the corresponding ra- 
tios are shown; compare legend, (c) Simulated constant-Q-scan ('scan 
A'), (d) Ratio of the intensity obtained in the simulation of 'scan A' in 
configuration parab. #2b versus elliptic #1 , Straight and parab. #2. (e) 
Simulated constant-E-scan at 2 meV ('scan C) for different models, 
(f) A simulated w-scan at (1 0). 



tion. 

Finally, we note that a reduction of the Ac/-shift to- 
wards zero does not yield any improvement but even a 
clearly diminished intensity for 'scan C in Fig.[T2l(e). 

7. Conclusion 

We have performed Monte-Carlo simulations in or- 
der to determine the optimal focusing design for a cold- 
neutron triple-axis spectrometer which will be installed 
at the end of a curved neutron guide (NL1 at the FRM- 
II). The elliptic guide 13111 turns out not to be the ideal 
choice since the elliptic configuration focuses also neu- 
trons with higher divergence into its focal point. In con- 
trast, the parabolic alternative yields much higher inten- 
sities by focusing neutrons with small divergence into 
its focal point. 

The highly focusing design of the parabolic con- 
cept also results in a distinctly better energy resolu- 
tion and does not suffer that much from a severe multi- 
peak structure transverse to Q like the elliptic concept. 
Hence, intensity as well as energy and Q-resolution are 
superior in the parabolic concept. This improved per- 
formance of parabolic focusing could be also verified 
in simulations for phonon-scans and of scans across the 
elastic line of an incoherent scatterer. 

However, for steep dispersions and low neutron ener- 
gies, both focusing concepts either elliptic or parabolic 
become quite similar having a significant intensity dis- 
advantage compared to a standard straight guide. 

The monochromator is the essential element of the 
whole focusing concept and can be optimized for the 
required curvatures in the double-focusing operation 
mode. If all distances are small and the curvatures be- 
come large, a smaller crystal size and a slightly en- 
hanced crystal mosaic are beneficial in order to follow 
the large curvatures and probably also to transport a 
larger part of the divergence bandwidth provided by the 
guide. 

We may conclude that for an instrument at the end of 
a neutron guide, the parabolic concept is clearly supe- 
rior to the elliptic concept. In most cases the parabolic 
concept is distinctly better than a standard straight guide 
unless the measured dispersion is steep and high Q- 
resolution is desired. 
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